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Abstract 

A wide field of view Cherenkov/fluorescence telescope array is one of the main components of the Large High 
Altitude Air Shower Observatory project. To serve as Cherenkov and fluorescence detectors, a flexible and mobile 
design is adopted for easy reconfiguring of the telescope array. Two prototype telescopes have been constructed and 
successfully run at the site of the ARGO-YBJ experiment in Tibet. The features and performance of the telescopes 
are presented. 
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1. Introduction 

The energy spectrum of primary cosmic rays spans 
almost 12 orders of magnitude, from 10^ eV to 10^' eV, 
and can be well fitted by a simple power law except 
in several small energy regions. A region called the 
"knee" of the spectrum existing at around 10'^ eV is one 
of these regions where the spectrum becomes steeper 
at higher energy side. Many experiments have ob- 
served this phenomenon; however, controversial argu- 
ments on its origin persist because of limited discrim- 
ination power on the primary cosmic ray composition 
and ambiguities in nucleus-nucleus interaction model- 
ing. These two aspects are closely related to each other 
Modern balloon borne experiments, such as ATIC |1] 
and CREAM |2|. have efficiently measured the energy 
spectra of individual elements at the top of the atmo- 
sphere. The energy spectra for all nuclei are measured 
up to ~100 TeV which is not far from the "knee". 
Because the detector area is constrained by the pay- 
load, the spectrum measurement has to be extended to 
a higher energy using a ground based air shower detec- 
tor array. The spectrum should initially be measured 
well below 100 TeV to create an overlap with the bal- 
loon experiments which serve as absolute calibrations 
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for the ground-based techniques. Identifying the in- 
dividual components of cosmic rays continues to be a 
major challenge in ground-based experiments. Mul- 
tiple parameter measurements on an air shower seem 
to be a plausible approach. The ultimate goal is to 
separate individual species out of the total observed- 
event samples and measure a clear individual "knee" 
for every single species, enabling the discovery of the 
origin of the "knee". As one of the major scientific 
goals of the Large High Altitude Air Shower Obser- 
vatory (LHAASO) project [jsi&l], the energy spectrum 
for a separated composition will be measured at ener- 
gies above dozens of TeV. To tag each primary par- 
ticle that causes an air shower, the atmospheric depth 
of the shower maximum should be measured as one 
of the important parameters. The wide field of view 
Cherenkov/fluorescence telescope array (WFCTA), one 
of the components of the LHAASO project, is designed 
to accomplish this goal. 

A portable design of WFCTA telescopes is adopted 
to maximize the flexibility of changing the configura- 
tion of the array of telescopes. The elevations, point- 
ing directions, and locations of the telescopes are then 
easily reconfigured. This is one way of using the same 
telescopes to serve as both fluorescence and Cherenkov 
detectors. In the fluorescence detector, the telescopes 
are tilted down to a horizontal position. In such an 
operational mode, which is analogous to the HiRes 
experiment OlSll . most of the Cherenkov photons are 
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avoided except those that are scattered onto the field 
of views (FOVs) of the telescope, such as in the fluo- 
rescence detector of the telescope array experiment 1 10] 
and the fluorescence detector of the Pierre Auger experi- 
ment IMS]- Only the fluorescence light from the shower 
is collected together with the scattered Cherenkov light 
to trigger the telescopes. This requires showers hav- 
ing much higher energy, usually above 100 PeV, such 
as in the HiRes prototype experiment [ISJ . because 
the fluorescence light by a single electron is consid- 
erably weaker and isotropic. In the Cherenkov detec- 
tor, the telescopes run in high elevation mode to di- 
rectly measure Cherenkov light from the showers, sim- 
ilar to what was done in the Dice experiment [16]. A 
Cherenkov light radiation provides considerably more 
photons along the shower axis that are useful for lower- 
ing the shower energy. 

In 2007, two prototype Cherenkov telescopes jsl 0] 
were deployed at Yangbaijing (YBJ) Cosmic Ray Ob- 
servatory near the ARGO-YBJ experiment |7]. More- 
over, two WFCTA telescopes have been successfully 
running in Cherenkov mode beginning August 2008. To 
date, millions of cosmic ray events that simultaneously 
trigger the telescopes and the ARGO-YBJ detector car- 
pet array have been collected. An analysis of these 
events is carried out to study the performance of the tele- 
scopes. Detailed descriptions of the telescopes and the 
analysis of the findings are presented in this paper. 

Several details about the apparatus are presented in 
Section 2. The detector calibration is then discussed in 
Section 3. The test-run of the two telescopes and re- 
sults are reported in Section 4 including summaries on 
the detector performance. The conclusions drawn are 
provided in the last section. 

2. Apparatus 

The two prototype telescopes are deployed near the 
ARGO-YBJ carpet detector array at a longitude of 
90.53°E, and a latitude of 30.11°N and 4300 m a.s.l. 
One telescope is about 25 m away from the west side 
of the ARGO-YBJ array. The other is also 25 m away 
from the south side of the array with separation distance 
between the two telescopes is 50 m. Each telescope has 
an FOV of 14° in elevation by 16° in azimuth. The fo- 
cal plane camera is made of a 16x16 photomultiplier 
tube (PMT) array, and the pixel size is approximately 
1°. Because both telescopes are tilted up to 60° pointing 
in the same direction, they can be operated in stereo- 
scopic mode, i.e., showers striking an area covered by 
the telescopes will be seen simultaneously. Since the 
Cherenkov light from a shower is very concentrated in a 



104 forward region; thus, the telescopes can be triggered by 

105 showers coming within a cone of approximately 8° with 

106 respect to the main optic axes of the telescopes. 

107 The entire telescope system is composed of an op- 

108 tic ultraviolet light collector, a focal plane camera, front 

109 end electronics (FEE) based on 50-MHz flash analog- 
10 to-digital-converters (FADC), data acquisition (DAQ) 
Ml based on an embedded ARM processor and PC 104 bus, 
12 power supplies for low and high voltages, and a slow 
113 control system. Everything is installed in a shipping 

14 container with dimensions of 2.5 mx2.3 mx3 m (FiglTJ. 

15 Mirrors are mounted at one end of the container and 

16 the camera is located at the other end where the focal 

17 planes of the mirrors are. The FEE and DAQ are placed 

18 at the back plane of the PMT camera. A glass window 

19 is installed at the entrance aperture to keep dust from 

120 entering the apparatus. The container is mounted on a 

121 dump-truck frame with a hydraulic lift that allows the 

122 container to be lifted up from 0° to 60°. The mobility 

123 of the entire telescope allows for freely switch between 

124 configurations of the telescope array for different obser- 

125 vational modes. The architecture of the electronic data 

126 acquisition and the slow control system are shown in 

127 Fig|2] whereas that of a sub-cluster is shown in Fig[3] 

128 The PMT signals are processed using an analog process- 

129 ing board (AB) and then a digitization board (DB). The 

130 first level trigger (FLT) is generated in the DB on a sub- 

131 cluster. After the FLT is determined, 256 FLTs are then 

132 sent to the trigger board (TB). The second level trigger 

133 (SLT) and the third level trigger (TLT) are then deter- 

134 mined in the TB. After this, the event trigger from the 

135 TB is fanned out by a bus driver board (BDB) and sent 

136 back to each of the DB and GPS boards. The data are 

137 initially stored in the bufler of the DB when the DB has 

138 received the event trigger, after which the data that in- 

139 elude the GPS time are read using TS7200. Finally, the 

140 data are stored in a PC in the laboratory via Ethernet. 

141 A detailed description of the detector, divided into the 

142 following 8 subsystems, is provided: 1) optics, 2) cam- 

143 era, 3) FEE, 4) trigger system, 5) DAQ, 6) power supply 

144 system, 7) slow control system that includes monitoring 

145 of everything, and 8) calibration. 

146 2.1. Optics 

147 A 4.7 spherical aluminized mirror, composed of 

148 20 hexagon-shaped segments, is used as an ultraviolet 

149 light collector Each segment is subjected to a strict con- 

150 trol of the surface quality and their geometrical and op- 

151 tical properties. The reflectivity is greater than 82% for 

152 light having a wavelength >300 nm. The radius of the 

153 curvature of the segments is 4740 mm with a tolerance 

154 of ±20 mm. 




Figure 1 : Photograph of the telescope with the doors open. 
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photons can be efficiently determined by simply taking 
the average of the directions of the registered pixels as 
weighted by the measured charge in a pixel. Contrast- 
ingly, the size of a spot is optimized to be similar to 
pixel size to avoid sharing of incident photons by an ex- 
cessive number of pixels. 

Because of the abbreviation of the spherical reflector, 
the spot size changes across the focal plane in a rather 
large FOV. Large comas also occur at large off-axis an- 
gles. The uniformity of the spot size over the camera 
is also optimized by locating the camera slightly away 
from the focus. A distance of 2305 mm between the 
mirror and the camera is eventually set according to a 
detailed ray-tracing calculation. 

Each mirror segment is mounted on a spherical steel 
frame with three adjustable screws. The pointing ori- 
entations of all segments are adjusted toward the geo- 
metric center of the curvature. Using a laser beam, the 
orientations are calibrated to be less than 7.6 arc seconds 
from the nominal direction. 



50 MHz 
clocks 



Sub-cluster 



Sub-cluster 



16 Sub-clusters 



256 FLTs 



Event 
triggers 



Bus 



SLT && TLT 
Generator 



Clock 
Distribution 



! Event 
i trigger i 



PC 
In Lab 



TB 

Event trigger 
GPS Board 



Bus Driver T 
Trigger Driver 
& Distribution 



MOOMbps I-ioOm,„, 
Ethernet « 



Bus 



TS7200 
(PC104) 



193 
194 




Figure 2: Communications diagram of one telescope; for details of 
the sub-cluster see Fig|3] 



The size of a light spot on the focal plane where the 
PMT camera is located is designed to be similar to the 
size of a pixel. The sensitivity of a PMT is non-uniform 
across the surface of the photocathode |17], and small 
gaps exist between PMTs in the camera. For a spot 
that is exceedingly larger than the dimension of the gaps 
and the typical spatial scale of the non-uniformity across 
the photocathode, light from a specific direction will be 
shared by a few adjacent PMTs. Thus, effects stemming 
from the overall non-uniformity of response across the 
entire camera are reduced. The direction of incident 
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2.2. PMT Camera 

Photons in a spot at the focal plane are recorded by a 
camera composed of 256 pixels which are 40 mm Pho- 
tonis hexagonal PMTs (XP3062/FL). The camera is en- 
closed in a box and mounted on a frame at a distance of 
2305 mm from the spherical mirror. 

The PMTs in a telescope are operated at a gain of 
6x10^. To achieve the best uniformity at hardware level, 
a resistor is placed between the high voltage power sup- 
ply and each PMT base to compensate for the gain dif- 
ference between PMTs. The gains (G) for all PMTs and 
their responses to the supplied voltage, i.e., G oc y^^, 
are calibrated. All yS values are measured and recorded 
in a database in the laboratory for future use. The PMTs 
are then sorted according to their gains and grouped into 
two classes. For instance, the working voltages of the 
PMTs of one of the telescopes distribute 1088 V to the 
telescope in a common power supply with a variance 
of 73.97 V. The average fi of these tubes is 5.9 with a 
variance of 0.6. 

The 16 PMTs, as an integrated unit, are soldered on 
a high voltage board (HVB) that distributes a negative 
high voltage to all cathodes and dynodes of the PMTs. 
The voltage division scheme is recommended by Pho- 
tonis to yield the maximum gain of the tubes. The max- 
imum output current of the PMTs allows a range of 3.5 
orders of magnitude in which the non-linearity of all 
tubes is less than 8% according to a calibration in the 
laboratory iITtIi. The anode signals are finally DC cou- 
pled to the FEE. 
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Figure 3: Photograph of a sub-cluster (left) and schematic of the sub-cluster (right). 
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2.3. FEE and Digitization 

The FEE is located behind the HVB, which is com- 
posed of an AB and two DBs, to avoid a long distance 
transmission of the analog signals. Such a module is re- 
ferred to as a sub-cluster (see the left figure of Fig|3]l in 
this paper. A block diagram on the procession of sig- 
nals from the PMTs to the on-board data storage chip is 
shown in the right figure of Figl3] The signal of each 
PMT is transmitted to the first amplifier on the AB for 
noise filtering and pulse stretching. Upon division into 
high or low gain channels, the signal is then sent to the 
DB for digitization and further processing. The rest of 
this section describes the processing procedure in detail. 

2.3.1. Analog Processing Board 

Each AB has 16 channels for the signals from 16 
PMTs on HVB alone. Each channel has a shaping cir- 
cuit and amplifiers for low or high gains. The board has 
four main functions, namely: 

• low-pass filter of 20 MHz, 

• expanding narrow pulses, 

• dual-gain system for covering a dynamic range 
over 3.5 orders of magnitude, 

• and receiving PMT signals from HVB and per- 
forming a single-ended-to-difFerential conversion. 

The anode signals are fed to a four-pole low-pass 
filter based on an AD8039 through DC coupling. 
Cherenkov photons generated by all shower particles 
move at almost the same speed as the charged particles; 
thus, all photons generated over the entire shower de- 
velopment beginning from the top of the atmosphere hit 
the cathodes at almost the same time. The duration of 
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Figure 4: Stretching ratio as a function of the input pulse duration. 



the pulse lasts only a few nanoseconds, which is shorter 
than the typical response time of the PMT, e.g., a pulse 
duration generated by a single-photo-electron (SPE) is 
approximately 12 ns on average. Therefore, such a nar- 
row pulse has to be stretched (e.g., at > 50 ns) to be 
able to measure the charge at a sampling rate of 50 
MHz, as preselected for the FADCs on DBs. Taking 
into account the optimized bandwidth of the noise filter, 
the stretching ratio is selected in such a manner that the 
narrower pulses are stretched further, and pulses wider 
than 120 ns are essentially not stretched. In FigH) the 
stretching ratio is plotted as a function of the pulse du- 
ration. According to this, the original waveform is re- 
constructible, and the timing of the pulse can be cor- 
rected with an uncertainty less than the duration of the 
stretched pulses. 

The amplification of the PMT signals with a gain of 
2.67 is another feature of the shaping circuit. The value 
of the gain is also optimized together with the filter and 
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Figure 5: A Cherenkov signal with sky background; the first 100 bins 
are low gain background and the remaining 800 bins are high gain 
background and the Cherenkov signal. The Cherenkov signal is in the 
current frame (from the 300th bin to the 600th bin). 

266 the stretcher. The AD8039 has a sufficiently high speed 

267 (350 MHz), low power dissipation, low cost, low noise, 

268 low distortion, and a nonlinearity less than 1 % in a range 

269 of the input from 0.5 mV to 800 mV, which nearly fits 

270 the entire range of the PMT signals. 

271 To maintain good linearity over a wide dynamic range 

272 of 3.5 orders of magnitude in charge, a dual gain system 

273 is designed. Signals coming out of the AD8039 are split 

274 into two channels and are then separately amplified by 

275 two AD8138. With such a design, considerable flex- 

276 ibility in obtaining different resolutions of the charge 

277 measurements for pulses with different pulse heights is 

278 achievable by choosing the range covered by the high 

279 gain channel. In this paper, a ratio of the gains of 1:8 is 

280 selected as a result of the optimization between the dy- 

281 namic range and the resolution in the charge measure- 

282 ment. The nonlinearity of AD8138 is less than 2% at 

283 both gains of 1 and 8 in the entire range of the input 

284 signals. 

285 Another excellent feature of the AD8138, the con- 

286 version between a single-ended input and differential 

287 output for instance, makes it the best choice for an 

288 amplifier. Superposing an offset to the output of the 

289 AD8138 as a pedestal of the signal before feeding into 

290 the positively-polarized FADC is also a convenient ap- 

291 proach. To manage possible undershoot of the pulses 

292 (see Fig|5ll, a pedestal is preset higher than the under- 

293 shoot. 

294 2.3.2. Digitization Board 

295 Each DB has 8 pixel channels, which includes 16 

296 FADC modules managed by two field programmable 

297 gate array (FPGA) modules. The 10-bit FADCs, 

298 AD9215, has a good cost-performance ratio with a lin- 



299 earity better than 2% over a range of 40 to 900 counts. 

300 However, it is not sufficient to cover the required dy- 

301 namic range of the detector Two FADCs are then used 

302 for the dual gain system. At an input range of 2 V, the 

303 FADC continuously digitizes analog signals with a sam- 

304 pling rate of 50 MHz (20 ns/bin) and a resolution of 2 

305 mV per count. In the dual gain system, a range of to 

306 500 photoelectrons with a resolution of 1 .7 counts per 

307 photoelectron is set for the high gain channel, whereas 

308 a range of to 4000 photoelectrons with a resolution 

309 of 0.21 counts per photoelectron is set for the low gain 

310 channel. 

311 A digitized waveform is collected by an FPGA, Xil- 

312 inx XC3S1000, and fed into a pipeline with a length 

313 of 1500 clock cycles. Every 300 cycles is defined as a 

314 frame in which a single channel trigger is formed. Such 

315 a long pipeline enables an enduring waveform wait- 

316 ing for global trigger formation and transmission. This 

317 pipeline also allows three frames, i.e., previous, current, 

318 and post ones, to be recorded once a global trigger is 

319 received. This guarantees storage of a complete wave 

320 form regardless of when the pulse starts in the "current" 

321 frame. 

322 The FPGA also makes a choice between the signals 

323 from high gain and low gain channels when the bit 

324 stream flows in. Any signal higher than 900 counts trig- 

325 gers a switch from the high gain channel to the low gain 

326 channel. This way, an effective overall dynamic range 

327 of 12 to 13 bits is achieved. 

328 2.4. Trigger system 

329 A final trigger among the telescopes is determined 

330 based on the three trigger levels, namely, the single 

331 channel trigger (the lowest level), telescope trigger (the 

332 second level), and event trigger (the highest level). The 

333 topology of the three-level trigger algorithm is outlined 

334 in Fig|6] with corresponding details discussed in follow- 

335 ing subsections. 

336 2.4.7. First Level Trigger in a Single Channel 

337 The FLT is formed in the FPGA located on the 

338 DBs provided that the signal-to-noise ratio in a win- 

339 dow is greater than a given threshold (e.g., 4 as a typi- 

340 cal value). The width of the window is predetermined 

341 for corresponding observation modes (e.g., 8 bins for 

342 the Cherenkov light signals). Running over the entire 

343 frame of 300 cycles once in a bin, 293 sums of FADC 

344 counts in the windows (denoted as WINSUM) are pro- 

345 duced. The average of WINSUM is calculated using 

346 fliese WINSUMs avoiding the maximal WINSUM and 
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Figure 6: Three-level trigger system. SM is stereo or mono enable 
signal, GL is global or local enable signal, Tl is the telescope trigger 
from telescope No. 1 and T2 is from telescope No. 2. See detailed 
explanation of the diagram in the text. 

347 18 WINSUMs on each side of the maximum. The stan- 

348 dard deviation (cr) of the WINSUM around the aver- 

349 age serves as a measure of the background fluctuation. 

350 If the maximum of WINSUMs exceeds the average by 

351 ncr, the FLT is formed, where n is preselected from a 

352 Ust of VT2, VT6, V24, V32 and V40 before starting 

353 the run. All trigger signals from 64 FPGAs are encoded 

354 into bitmaps and transmitted to the Trigger Board (TB) 

355 in parallel (see Fig|6l). 

356 The frame size of 300 bins is selected for a full scale 

357 array with more than 24 telescopes covering an area of 

358 1 km^, over which a highly inclined air shower takes 

359 several microseconds to cross in the fluorescence light 

360 observational mode. Such a frame is large enough to 

361 fully contain the entire shower; it is also large enough 

362 for Cherenkov light. All Cherenkov photons arrive at 

363 the telescopes at almost the same time (few ns). There- 

364 fore, an air shower signal appears only within 5 or 6 

365 bins after shaping and the night sky background is mea- 

366 sured in the rest of the frame (see Fig|5]l. Such contin- 

367 uous measuring of the sky background is highly useful 

368 not only in estimating the signal-to-noise ratio but also 

369 in monitoring the transparency of the atmosphere using 

370 well-known bright stars in the field of view. 

371 2.4.2. Second level trigger for a single telescope 

372 A trigger for a telescope registered by an air shower 

373 (denoted as second level trigger or SLT) is formed when 

374 a specific pattern of triggered PMTs corresponding to 

375 a possible air shower is found in the camera of a tele- 

376 scope. A pattern recognition technique, developed in 

377 the Pierre Auger experiment 1 11 ] that is operated in a 

378 FPGA located on the TB, is applied. Numerous pat- 

379 terns, such as a fully filled circle formed with one hexag- 

380 onal pixel surrounded by six others, or a straight line 

381 formed by six aligned pixels, are pre-loaded into the 



Figure 7: Two typical patterns in the second level trigger. 

382 FPGA as a look-up table. Once all 256 FLT signals are 

383 collected, the FPGA matches all the pre-stored patterns 

384 with the observed one within a 6 x 6 matrix of pixels. 

385 Moreover, it keeps such a box running throughout the 

386 entire PMT camera with a step of one row or one col- 

387 umn. The SLT is formed as long as any one of the pre- 

388 stored patterns is matched. According to the simulation 

389 for air showers, mainly two types of patterns exist. Flu- 

390 orescence light images of air showers seen from several 

391 kilometers away tend to form line-shaped patterns on 

392 the camera, whereas Cherenkov light images of showers 

393 hitting the telescope head-on tend to form round-shaped 

394 patterns, as shown in Fig|7] There are 16 round-shaped 

395 patterns and 729 line-shaped patterns in a 6 x 6 box. 

396 To speed up the formation of SLT, the pattern com- 

397 parisons in the box are conducted in parallel. All 121 

398 bitmaps of the boxes are generated by sliding the box 

399 are done in parallel. All 121 bitmaps of the boxes are 

400 generated by sliding the box throughout the camera. A 

401 pattern comparison algorithm is then carried out among 

402 the 121 bitmaps. The telescope trigger (SLT) is formed 

403 in 123 clock cycles, i.e., 2.46 fis. 

404 2.4.3. Third level trigger for entire array of telescopes 

405 An event trigger for the entire array of telescopes 

406 (denoted as the third level or TLT) is generated using 

407 one of the FPGAs used for SLT. For the two proto- 

408 type telescopes, only two modes, namely, stereoscopic 

409 and monocular observation of showers, exists if they are 

410 configured in such a way that the two telescopes have a 

411 maximal overlap of the FOV. In stereoscopic mode, the 

412 two telescopes are required to simultaneously observe a 

413 shower. In monocular mode, each of the two telescopes 

414 can trigger the entire site. The operational mode should 

415 be selected at the beginning of a run by assigning two 

416 controlling parameters, SM and GL, as marked in Fi^ 

417 The two telescopes communicate with each other 

418 through a two 80 m coaxial cables by sending telescope 



6 



419 
420 
421 
422 
423 

424 
425 
426 
427 
428 
429 
430 
431 
432 
433 
434 
435 
436 
437 
438 
439 
440 
441 
442 
443 
444 
445 
446 
447 
448 
449 
450 
451 
452 
453 
454 

455 
456 
457 
458 
459 
460 
461 
462 
463 
464 
465 
466 
467 
468 



triggers (SLT) out and receiving the event trigger (TLX). 
Such a trigger function is also useful for receiving an 
external trigger and delivering a trigger to other co-sited 
detectors. This makes the telescope an open and modu- 
larized system. 

2.5. DAQ system 

The philosophy of the design of the Cherenkov tele- 
scope DAQ is integrated and compact so that the en- 
tire DAQ can be arranged on the backboard of the focal 
plane camera. It allows for the maximal mobility of the 
telescope for flexibility in switching between configura- 
tions. Because of the limited space, a low power con- 
suming and compact-embedded online computer based 
on an ARM processor, industrial standard PC 104 bus, 
and flash disk are selected as the backbone of the DAQ 
system, eliminating the need for moving parts such as 
the CPU fan and hard drives. All of the components are 
integrated on a bus driver board, which bridges the data 
storage disk and the 32 DBs that are connected by flat 
cables. Because the event rate is not extremely high, 
i.e., around 1 Hz, a band width of 14 MHz for a PC104 
bus is sufficient for transferring every bit of the 6 //s long 
wave forms to the disk for all 256 channels. It is also ex- 
tremely advantageous for system debugging. Each tele- 
scope is an independent detector with a complete DAQ 
system on board. The communication with the rest of 
the experiment is through a 10/100 Mbps Ethernet. 

Each telescope has its own independent DAQ, devel- 
oped using C-n- language under Linux. Operating in 
a polling mode instead of an interrupt mode, the DAQ 
checks the interface status regularly for whether the data 
are ready or not. If so, the data in the buffer are read 
and immediately stored in the disk. Before a night shift 
ends, the data in the hard disk are moved to the Institute 
of High Energy Physics (IHEP) in Beijing for further 
analysis. 

2.6. Power supply, slow control, and status monitoring 
The two power supplies with +7 V (maximum current 

80 A) and -7 V (maximum current 20 A) are installed in 
the container. Regulators are used at each AB and DB 
to further stabilize the voltage. An adjustable HV power 
supply with a maximal output of -2000 V and 100 mA 
is used for each telescope. 

One of the difficulties at high altitudes is the heat dis- 
sipation of the power supplies and electronics enclosed 
in a metal box. The total power consumption is about 
50 A at +7 V and 16 A at -7 V. A forced cooling sys- 
tem is necessary in the prototype experiment to improve 
the heat dissipation of the entire camera and the power 
supplies. 
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The entire telescope system is powered by an unin- 
terrupted power supply with a sufficient battery backup. 
It protects the telescopes from damage when a blackout 
occurs. 

The detector is designed to work in remote control 
mode, including the opening and closing of the doors, 
turning the power supplies on/off, enabling/disabling 
the high voltage (HV) power supply and low voltage 
(LV) power supplies, and switching the UV LEDs on/off 
for detector calibration. All kinds of controls are real- 
ized using the on-board computer through a COM port. 
Carefully monitoring the status of each telescope is nec- 
essary, including the door status, voltages, and tem- 
peratures at different places (e.g., enclosed areas inside 
the camera, backboard of the PMT camera, and inside 
the housing of the UV LED). All parameters are mea- 
sured and recorded through an 8-channel 12-bit analog- 
to-digital-converter on the board of the embedded com- 
puter, TS7200. 

All of the controls are performed through a user in- 
terface running at the embedded computer connected 
through Ethernet from an operational center, 3000 km 
away in IHEP in Beijing. A small portion of the data 
can be copied to display the event during the operation. 

3. Calibration 

3.1. Method of calibration 

An accurate shower reconstruction requires a con- 
verting factor from a pulse area in terms of FADC 
counts to the number of photons for each pixel. There- 
fore, the absolute calibration of the detector response is 
essential. To achieve this objective, having an accurate 
knowledge about PMT cathode effective areas, cath- 
ode quantum efficiency, PMT gains, amplifier gains, 
and digital converting factors, is necessary. Measur- 
ing all these effects item by item is difficult. In this 
paper, a method similar to the HiRes experiment lfl4ll 
and Pierre Auger experiment lll2ll calibration procedure, 
which considers the entire effect, is applied to measure 
the overall response of each pixel. The procedure is dis- 
cussed in the paragraphs that follow. 

Being mounted at the center of the mirror, a UV-LED 
(375 nm) light source with a diffuser is used for the cal- 
ibration of PMTs in the camera by beaming nearly uni- 
form light to every pixel. The LED light density on the 
camera surface is calibrated using a pre-calibrated probe 
detector located beside the camera. The calibration of 
all pixels in the cameras is performed twice a day, i.e., 
before and after the daily operation. The crucial part is 
the measurement of the absolute number of UV photons 
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at each PMT cathode from a pulse emitted by the LED, 
which is performed in two steps. 

First, we move one PMT over all the places on the 
frame of the PMT camera to measure the uniformity of 
the LED light density on the camera surface. This has 
to be accomplished before the PMT camera is installed. 
The light density is a function of the polar angle, 9, in a 
form of cos'^0, where 6 is the angle between the connec- 
tion from the LED to the PMT location and the perpen- 
dicular connection from the LED to the PMT camera 
surface. The variation from the center to the corner of 
the PMT camera is within 7%. 

In the second step, the light density is measured using 
a pre-calibrated probe consisting of two XP3062 PMTs 
with the same FEE and DAQ as the two telescopes. The 
only difference is that the two PMTs in the probe are 
operated at a very high gain so that the single photoelec- 
tron can be measured. Therefore, the gain of the probe, 
Qprohe^ can be calibrated at any time. Then, the absolute 
calibration of the probe is performed at the HiRes lab 
at the University of Utah, USA, by comparing the re- 
sponse of a hybrid photo diode (HPD) pre-calibrated at 
NIST flS*] to the same light source. The probe and HPD 
are located side by side in front of a UV-LED at 355 nm. 
Using the HPD, we measure light density lu (number of 
photons per square millimeter) from the LED. We also 
measure pulse area F''^"'"' using the probe simultane- 
ously. At the operational site in Tibet, the light density 
from the LED mounted on the center of the mirror is 



calibrated to be Ij - II 



pprohe ^prohe 



/^camera 
^375 



j^camera 



COS 



7, where U and T stand 



for Utah and Tibet, respectively. 

Applying the knowledge obtained from the two steps, 
we have determined both the light density in front of 
the cathode and the pulse area in terms of the FADC 
counts for each PMT in the camera. Finally, the over- 
all converting factor for a pixel between the number of 
photons reached to the surface of the camera and cor- 
responding pulse area measured by the FADC counts 
behind the pixel is 



(1) 

hApMT cos^ff""""''' 
where Apmt is the geometric area of the PMT cathode, 
and C™"*"™ is the calibration constant for a pixel. The 
subscript 375 indicates that the calibration is done using 
UV light at 375 nm. The unit of the calibration constant 
is FADC counts per photon. 0'"'°'"^ and ff^"""^''" represent 
the angular locations of the probe and the pixel in the 
camera, respectively. A further correction according to 
the wavelength dependence of the quantum efficiency 
of PMTs is applied in the operation for the cosmic ray 
observation. 




54820 54840 54860 54880 54900 54920 54940 54960 54980 
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Figure 8: Absolute gain of one telescope from Dec. 2008 to May 
2009, as a function of time. 



3.2. Result of the calibration 

The probe was calibrated at the HiRes lab. The cali- 
bration results are shown in tablelT] the gain of the probe 
is 85.7+0.3 FADC count per SPE. The second column 
shows the lower LED photon density and the third col- 
umn shows the higher LED photon density that coincide 
with one another. 



Table 1 : The calibration results of the probe. 



/(7(photons/mm'^) 


0.359 ±4.8% 


0.538 ±4.8% 


^'^'''(FADC count) 


8959 +4 


13650 ±5 



The corresponding calibration constant for one of the 
telescopes is presented in FiglS] An average of the cal- 
ibration constants of all pixels in the camera is plotted 
with respect to time from December 2008 to May 2009. 
As mentioned above, the calibration constant is mon- 
itored every observational day. The systematic uncer- 
tainty of the calibration constant is estimated to be 7%. 
The downtrend of the constant shown in the figure indi- 
cates obvious decreases of gains in all newly produced 
of PMTs. The transmission of the glass window and re- 
flectivity of the mirrors are not take into account in the 
above calibration. These two effects will be monitored 
in so called end-to-end calibration using nitrogen laser 
in future. 

4. Performances and Results from Test Run 

4.1. Test Run Information 

The two telescopes began recording cosmic ray data 
in August 2008. Furthermore, both monocular and 
stereoscopic modes have been tested. About 500,000 



8 



594 coincidental events with the ARGO-YBJ experiment 

595 in stereoscopic mode and 700,000 coincidence events 

596 in monocular mode have been collected up to January 

597 2010. The average trigger rate is about 0.5 Hz in stereo- 

598 scopic mode and 0.7 Hz in monocular mode. 

599 4.2. Off-line Coincidence with the ARGO-YBJ Experi- 

600 ment 

601 All WFCTA telescopes and the ARGO-YBJ experi- 

602 ment recorded the cosmic ray arrival time based on a 

603 GPS. A time window of 8 //s containing a Cherenkov 

604 event is searched for coincidence with the ARGO-YBJ 

605 event stream, which is about 4 kHz. A difference be- 

606 tween the recorded event time by the two experiments 

607 for a matched event is typically less than 100 ns (Fig|9l). 

608 For a coincidence event, the shower geometry is mea- 

609 sured using the ARGO-YBJ detector. The distribution 

610 of shower arrival directions is shown in a 2-dimensional 

611 map, i.e., zenith angles versus azimuth angle, as in 

612 FiglTOl Approximately 85% of coincident events occur 

613 inside the FOV of the Cherenkov telescopes marked by 

614 the trapezium. The rest of the 15% of events have their 

615 images partially seen by the telescopes and sufficiently 

616 trigger the telescopes. The shower core distribution is 

617 shown in Fig[TT| The ARGO-YBJ experiment array and 

618 the two Cherenkov telescopes are marked as the rectan- 

619 gle and two dots in the figure, respectively. For events 

620 that have the reconstructed cores inside the ARGO-YBJ 

621 array, the shower parameters, such as the number of hits 

622 on the carpet detector (A^/,„) and shower geometry, are 

623 well measured. A distribution of Nhti of these events is 

624 plotted in Fig [12] According to the number of hits as a 

625 function of the primary energy of proton from ARGO- 

626 YBJ detector 119], the mode energy of protons is about 

627 40 TeV. This estimates the threshold of two Cherenkov 

628 telescopes in stereoscopic mode. 

629 In Fig [13] two Cherenkov images are shown in the 

630 event display. The top image is caused by an event at 

631 a 10 m Rp (the impact parameter of the shower to the 

632 telescope) and the bottom image stems from an event at 

633 a 173 m /?p. The image is clearly more elongated for the 

634 farther event. Using the well-defined image parameters 

635 created by Hillas [21], length and width, this effect is 

636 more quantitatively presented in Fig[T4l The widths of 

637 the images seem to be no longer shrinking once show- 

638 ers are sufficiently far from the telescopes (e.g., farther 

639 than 100 m). The showers essentially resemble linear 

640 patterns. 

641 4.3. Electronic noisy and sky background 

642 For a triggering system completely based on the 

643 signal-to-noise ratio, such as the WFCTA telescopes. 
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Figure 9: The difference time between WFCTA and ARGO in coinci- 
dence event. 




azimuth (deg.) 

Figure iO: Distribution of coincidence event.s are seen by ARGO and 
Cherenkov telescopes simultaneously over zenith and azimuth angle. 
About 85% of coincidence events are located in the Cherenkov tele- 
scope FOV (trapezium). The image of the rest of the 15% of coinci- 
dence events is partially in the Cherenkov telescope FOV. 

644 the sensitivity is constrained by the noise level. For a 

645 most optimized system, the electronic noise must be 

646 negligible compared with the night sky background. 

647 Both of them are measured during the operation as the 

648 door is closed and opened, respectively. 

649 The sources of electronic noise include thermal 

650 noises in the PMTs, noises from the HV power supply 

651 and LV power supplies, noises of amplifiers, and finally 

652 from the counting error of FADCs. The PMT Photonis 

653 XP3062/FL has a very low dark current; therefore, the 

654 thermal noise level is sufficiently low so that it can be 

655 ignored compared with the other sources of electronic 

656 noises. The HV power supply has a ripple less than 

657 0.02% in terms of RMS of rated voltage. This con- 

658 tributes a variation within 0.12% in terms of the gain 

659 of the PMTs. Voltage regulators are used to block pos- 




x(m) 



Figure 11: Distribution of shower core. Two Cherenkov telescopes 
(two dots) and the ARGO-YBJ experiment array (rectangle) are also 
marked in the figure. 








Iog10(number of hits) 

Figure 12: Distribution of ARGO nhit of these coincidence events, 
whose reconstructed core is located in the ARGO-YBJ cluster array. 

660 sible noises from the LV power supplies. All ampli- 

661 fiers, ADS 138 and AD8039, are selected to generate 

662 a very low noise level. The noise from the FADC is 

663 about 0.5 FADC counts per tube. Taking into account 

664 all additional sources such as the distributing capacity 

665 on boards and connectors, the total electronic noise, in- 

666 eluding PMT and high voltage power supply, is typi- 

667 cally about 1.0 FADC counts per tube. It is measured 

668 during the calibration with the LED (FigfTSll as long as 

669 the signal is avoided. 

670 The night sky background is measured similarly as 

671 a shower trigger is formed when the doors are open 

672 (FiglSj). The fluctuation in the night sky back-ground is 

673 typically more than 2.2 photon electrons per 20 ns per 

674 tube on a clear moonless night. This suggests a much 




Figure 13: Two Cherenkov events; the Rp of the tope event is 10 m 
and that of the bottom event is 173 m; the image of the bottom event 
has a longer tail than the top event. 

675 Stronger noise than electronic noises. One of the impor- 

676 tant sources of the sky background is light from the stars 

677 (FigfTSll. which shows the sky background at night, in 

678 which each pulse denotes stars passing though the FOV. 

679 A bright individual star can be traced when it passes 

680 through the telescopes. Stars provide numerous stable 

681 point-like sources that typically sweep across a tube in 

682 about 4 minutes. Those ideal point-like sources can be 

683 used for multiple calibration purposes. For instance, 

684 well-known bright stars can be used as light houses to 

685 establish the pointing direction of the telescope itself, 

686 with accuracy proved to be better than 0.05° ll20ll . A 

687 bright star can be seen as a very stable point-like source 

688 at infinity, which is a perfect tool for measuring the spot 

689 size produced by the optical system. The background 

690 of a tube signal steadily increases when the light spot 
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Figure 14: Length (filled circles) and width (open circles) as a func- 
tion of Rp. 



Figure 16: Sky background is monitored by a PMT in a clear moon- 
less night; each pulse corresponds to a star passing through the FOV 
of the telescope. 
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Figure 15: LED calibration for door closed. The low noise level of 
1 LSB of the complete electronics is to be seen before and after the 
LED pulse. 



moves into the FOV of the tube, and steadily decreases 
when it moves away from the FOV, forming a light pro- 
file. The light profile can be fitted using a Gaussian 
function, whose sigma denotes the spot size. The spot 
size grows larger when the star moves away from the 
center of camera (FigfTTb. 



5. Summary 

The telescopes were successfully run at YBJ from 
August 2008 up to July 2009. Millions of coincidence 
events with the ARGO-YBJ experiment have been col- 
lected. The performance of the telescopes was studied 
using these events. The trigger rate is about 0.5 Hz in 
stereo mode. Moreover, the mode energy of the tele- 
scope is 40 TeV when a pure proton composition is as- 
sumed. 
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Figure 17: A star moves away from the center of the camera, and 
its light profile recorded by PMT is drawn in filled circles. The light 
profile described by a Gaussian function is drawn in a solid line. The 
sigma of Gaussian denotes that the spot size is also marked. 



The features of the two WFCTA prototype telescopes 
are summarized as follows: 



708 • a 4.7 m spherical mirror, 

709 • a 16x16 PMT aiTay covers an FOV of 14° x 16° 

710 with 1° pixels, 

• dual gain system for a dynamic range to 3.5 orders 
of magnitude, 

• DC coupling and moduUzed design for electronics, 

• three-level online trigger logic: single channel trig- 
ger based on S/N ratio, telescope trigger based on 
pattern recognition, and event trigger for stereo- 
scopic observation. 
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718 • maximized mobility and the telescope can be up- 

719 Ufted from 0° to 60° in elevation. 

720 The absolute gains of the telescopes are calibrated us- 

721 ing calibrated LEDs mounted at the centers of the mir- 

722 rors. The systematic uncertainty of the calibration con- 

723 stant is about 7%. The pixel gains are monitored on a 

724 daily basis. 
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